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ABSTRACT 

We present global VLBI observations of the first-excited state OH masers in the massive star-forming region 
Onsala 1 (ON 1). The 29 masers detected are nearly all from the 6035 MHz transition, and nearly all are iden- 
tifiable as Zeeman pair components. The 6030 and 6035 MHz masers are coincident with previously published 
positions of ground-state masers to within a few milliarcseconds, and the magnetic fields deduced from Zee- 
man splitting are comparable. The 6.0 GHz masers in ON 1 are always found in close spatial association with 
1665 MHz OH masers, in contrast to the situation in the massive star-forming region W3(OH), suggesting that 
extreme high density OH maser sites (excited-state masers with no accompanying ground-state maser, as seen 
in W3(OH)) are absent from ON 1. The large magnetic field strength among the northern, blueshifted masers 
is confirmed. The northern masers may trace an outflow or be associated with an exciting source separate from 
the other masers, or the relative velocities of the northern and southern masers may be indicative of expansion 
and rotation. High angular resolution observations of nonmasing material will be required to understand the 
complex maser distribution in ON 1 . 



Subject headings: masers — ISM: individual (ON 1) 
lines: ISM 



1. INTRODUCTION 

Onsala 1 (ON 1) is a kinematically complex site of massive 
star formation. In the 1.6 and 6.0 GHz transitions of hydroxyl 
(OH), masers are seen in two disjoint velocity ranges: -2 to 
+6 kms" 1 in the north of the source and +11 to +16 kms" 1 
projected atop and to the south of the H II region dArgon et al.l 
20001 iFish et alJl2005t iNammahachak et aUl200a !Fishll2007t 
Fish & Reid 2007). A 4765 MHz excited-state maser at 



+24.1 km s -1 was also reported by lGardner & Martin-Pintadol 
(119831) but never redetected. The highly excited 13 441 MHz 
transition shows similar velocity structure as the 1.6 an d 
6.0 GHz masers dBaudrv & Desmursll2002l: IFish et al.ll2005l) . 
although an interferometric map of the 13 GHz masers has 
not yet been publish ed. The systemic v elocity of the H II re- 
gion as measured by Zhe ng et al.l (Q~985) in HI 6a emission is 
approximately +5 kms -1 , which led the authors to conclude 
that the redshifted masers seen atop the H II region are under- 
going infall. However, recent proper motion measurements 
strongly suggest that the masers are traci ng slow (~ 5 km s -1 ) 
expansion in ON 1 (Fish & Reid 2007 ), similar to what is 
seen in W3(OH) dBloemhof et alJll992Tl Wright et"aTJl2004at 
Fish & Siouwe rman 2007])] 

The ground-state OH masers in ON 1 hav e been observed 
with connected-e l ement interferometry (lHoetal.1 119831: 
Argon et all 120001: INamm ahacha k et all 120061: iGreenetalJ 
2007T and with very long baseline interferometry (VLBI) 
(IFish et al.l r 2005t IFish & Reidl 120071) on multiple occasions. 
However, the only VLBI obser vations of the excited-state 
6.0 GHz masers were taken by Desmurs & Baudrvl (119981) 
with a three-station array. Due to the poor image fidelity in- 
herent with such a sparse array, they were only able to detect 
7 masers at 6035 MHz and 2 at 6030 MHz, ranging in flux 
density from 0.5 to 7.1 Jy. The large increase in the number 
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of telescopes with 6.0 GHz capability has since made high- 
fidelity imaging of northern 6.0 GHz maser sources possible. 

A prime motivation for observing ON 1 with a high sensi- 
tivity, high angular resolution, high spectral resolution VLBI 
array derives from the experience observing W3(OH), an- 
other nearby massive star forming region. In W3(OH), it was 
found that the 6030 and 6035 MHz OH masers traced some 
areas where no ground-state masers have been observed and 
highlighted portions of other areas where ground-state masers 
are abundant, allowing a greater understanding of the large- 
scale ( Galactic cloud) structure del ineated by molecules in the 
source (Fish & Siouwerman 2007). Of the massive star form- 
ing regions visible from the north, ON 1 is one of the brightest 
6.0 G Hz OH maser source s, with numerous masers detected 
in the Baudry e t al.1 d 19971) survey. Thus, in addition to be- 
ing an interesting source per se, ON 1 is also an interesting 
test case to determine to what degree OH masers in the first 
excited state can provide information not available from the 
ground-state masers alone. 

2. OBSERVATIONS AND DATA REDUCTION 

ON 1 was observed on 2008 Jun 15 with a global ar- 
ray consisting of two National Radio Astronomy Observa- 
tory (NRAO) and seven European VLBI Network (EVN) tele- 
scopes with A = 5 cm capability (experiment code GS029). 
The participating stations were Effelsberg, the Green Bank 
Telescope (GBT), a single Expanded Very Large Array 
(EVLA) telescope, Hartebeesthoek, Jodrell Bank (Mark 2), 
Medicina, Onsala (85 ft), Torun, and Westerbork (single tele- 
scope). 

Both the 6030.747 and 6035.092 MHz main-line transitions 
of the first excited state (^3/2,/ = 5 /2) of OH were observed 
in dual circular polarization mode. For each transition, a 
500 kHz bandwidth (24.8 km s -1 ) was centered at an LSR ve- 
locity of +8 km s" 1 in order to cover the full range of velocities 
of maser emission previously detected in ON 1 . At correlation 
time, the data were accumulated for 1 s intervals and channel- 
ized into 512 spectral channels, corresponding to a spacing of 
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0.049 kms -1 . Both the parallel-hand and cross-hand polar- 
ization correlations were obtained in order to reconstruct all 
Stokes quantities. 

Data were reduced using the Astronomical Image Process- 
ing System (AIPS). The source 3C345 was used for instru- 
mental delay and bandpass calibration, and J2 136+0041 was 
used for polarization calibration. The polarization position 
angle of J2 136+0041 was taken to be -162.°68 from the 
average of C-band entries in the VLA/VLBA Polarization 
Calibration DatabaseQ. Visibility weights were modified by 
taking their fourth roots in order to minimize the extent to 
which the overall calibration solutions would be dominated 
by the Effelsberg-GBT baseline (and baselines to these two 
highly sensitive telescopes). Since OH masers usually dis- 
play substantial circular polarization, the ON 1 data were self- 
calibrated using one circular polarization of a bright maser 
spot in each transition, and the calibration solutions were ap- 
plied to both polarizations. 

The source J2003+3034 was used to determine the abso- 
lute position of a bright reference maser feature (listed in Ta- 
ble [TJ at each of the two frequencies. Images were produced 
in all four Stokes parameters as well as both right and left cir- 
cular polarization (RCP and LCP, corresponding to RR and 
LL correlations). Maser spots were identified in the RR and 
LL cubes where emission was above 6 times the rms noise in 
each velocity channel. For each maser feature, spectra were 
obtained by integrating over a tight rectangular box encom- 
passing the feature emission. 

3. RESULTS 

We detected 29 masers, of which the vast majority (26) are 
6035 MHz masers (Table Q}. The distribution of 6.0 GHz 
masers, shown in Figure Q] confirms at high er resoluti on the 
general distribution seen at lower resolution (Fish 2007). The 
spectra of recovered emission (total flux density in spots iden- 
tified as masers) are shown in Figure [2] The spec tra are very 
good qualitative matches to those obtained by iGreen et all 
(2002b. 

Neither 6035 MHz masers at vlsr ~ +1 or +2 km s" 1 in the 
north of the source nor any low-velocity 6030 MHz masers 
at all are det ected, despite their appearance in previous stud - 
ies of ON 1 dBaudrv et alJH997t iFish et all200aiFish]l2007l) . 
These masers were weak (peak flux density 0.2 Jy, or approx- 
imately the flux density of the w eakest detected maser in this 
work) in the Ba udry et al.1 (119971) Effelsberg observations. At 
least one such maser was seen above a brightness level of 
0.7 Jy beam -1 in the EVLA observations of lFishl (120071) . sug- 
gesting variability. We note a possible RCP 6035 MHz feature 
at (- 1 64, + 1 049) mas, vlsr = +2 .6 1 km s" 1 in our data, but the 
spot is seen only in one channel and only at 5.7 a, which does 
not meet our detectability criteria. 

Nearly every 6030 and 6035 MHz maser feature is found 
with an accompanying feature in the opposite polarization 
at a slightly different velocity, presumably due to Zeeman 
splitting in the presence of a magnetic field. There are 12 
or 13 readily identifiable Zeeman pairs constituting 24 or 
26 maser spots. The high Zeeman pairing efficiency (over 
80% of the detected 6 GHz masers can be paired) in the 
6 GHz transitions is consistent with the small Zeeman split- 
ting coefficients at 6030 and 6035 MHz as compared with 
the ground-state transitions. Observations of W3(OH) show 

3 |http : / /www .via ■ nrao ■ edu/astro/calib/polar| 
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FIG. 1. — Map of 6.0 GHz masers in ON 1, Squares and circles denote 
6035 MHz and 6030 MHz masers, respectively. Maser lines identified in RCP 
are shown as open symbols, while LCP masers are shown as filled symbols. 
Numbers indicate magnetic fields in milligauss as derived from Zeeman split- 
ting, with letters denoting the a ssociated feature s in Table [T] Contours show 
8.4 GHz continuum emission lArgonetalfi OOO). Offsets are measured from 
20'' 10'"09?0768, +3 1 ° 3 V 34'.' 923 (J2000). 
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FIG. 2. — Total recovered maser emission in ON 1. RCP is shown in bold 
and LCP in normal weight. The same vertical scale is used for both panels. 
Letters denote features in TablefT] 



a similar trend for increased Zeeman pairing ef ficiency with 
decreased Zeeman splitting in v elocity units (I Wright et al.l 
2004bi IFish & Siouwermarl2007l) . 

We detect linear polarization at the 4 a level or greater in 
several maser spots. At 6030 MHz, spot Z shows a linear 
polarization fraction (y Q 2 + U 2 /I) of 17% with an electric 
vector polarization angle of -32° east of north in the channel 
of pea k emission, approximately consistent with IGreen etaT] 
(2007). At 6035 MHz, spot K is 22% linearly polarized with 
a position angle of 52° and spot B is 17% linearly polarized at 
—19°. There is evidence that the position angle of these spots 
may vary across the emission line in frequency. There is a 
ridge of weak emission in eight consecutive velocity channels 
passing through spots Mm and Nn. The linear polarization 
in spots Mm and Nn is not strong in any individual velocity 
channel, but after integrating over all eight velocity channels 
in which emission is detected in this region, linear polariza- 



VLBI Observations of 6 GHz OH Masers in ON 1 



3 



TABLE 1 

Detected Masers at 6030 and 6035 MHz 





VLSR 


R.A. 


Decl. 


Flux Density 


Av 


Gradient 


Gradient PA 


Zeeman 


Pol. 


(kms -1 ) 


Offset (mas) 


Offset (mas) 


(Jy) 


(kms" 1 ) 


(mas (km s ) ) 


(°) 


Pair 




6035 MHz 


R 


12.36 


-883.778 


707.787 


0.201 


0.25 






A 


L 


12.61 


-883.134 


707.356 


0.262 


0.25 






a 


R 


13.82 


-855.931 


698.986 


2.499 


0.31 


1.40 


132 


B 


L 


14.06 


-855.852 


698.853 


1.374 


0.34 


1.15 


155 


b 


R 


13.92 


-847.432 


691.329 


0.405 


0.17 






C 


R 


15.16 


-257.680 


100.194 


0.622 


0.17 


7.72 


-70 


D 


L 


15.42 


-258.373 


99.586 


0.475 


0.18 






d 


R 


15.18 


-250.056 


102.927 


1.276 


0.26 


1.48 


-83 


E 


L 


15.45 


-250.018 


103.112 


0.824 


0.26 


2.09 


-41 


e 


R 


13.98 


-227.783 


102.493 


0.195 


0.20 






F 


L 


14.28 


-227.264 


101.794 


0.241 


0.19 






f 


R 


14.78 


-202.760 


87.994 


0.229 


0.24 






G 


L 


15.07 


-203.000 


87.858 


0.218 


0.24 






g 


R 


5.51 


-128.465 


994.596 


2.229 


0.26 


1.16 


-6 


H 


L 


5.5 b 


-128.492 


995.516 


0.227 








(h) 


R 


-0.20 


-49.760 


1013.784 


2.241 


0.27 


1.09 


-54 


I 


L 


0.48 


-49.883 


1013.685 


1.702 


0.27 


0.89 


-12 


i 


R 


-0.63 


-47.799 


1011.481 


0.668 




2.90 


127 


J 


R 


14.50 


0.000 


0.000 


10.274 


0.21 


0.75 


156 


K a 


L 


14.79 


0.148 


-0.515 


0.320 








k 


R 


13.67 


33.361 


533.130 


1.981 


0.25 


12.52 


-68 


L 


L 


13.76 


33.598 


533.325 


0.919 


0.26 


8.60 


-72 


1 


R 


14.40 


148.685 


-0.352 


0.583 


0.20 






M 


L 


14.47 


147.241 


0.679 


0.446 


0.22 






m 


R 


14.43 


160.529 


-6.392 


1.147 


0.21 






N 


L 


14.50 


160.754 


-6.787 


0.976 


0.22 


4.81 


130 


n 


6030 MHz 


R 


13.79 


-855.753 


699.107 


3.244 


0.28 


1.11 


169 


Z a 


L 


14.12 


-856.052 


698.982 


1.297 




1.12 


-124 


z 


L 


13.81 


-855.762 


699.078 


2.323 


0.18 






(z) 



NOTE. — Positions are offsets as measured from 20 h 10 n W0768,+31°3l'34"923 (J2000). 
a Reference feature. 
Velocity cannot be determined to an additional significant digit. See text. 



which it is detected. The positional gradients, given in Table 
Q] indicate the magnitude and position angle of the gradient 
(increasing in the direction of positive velocity). Maser com- 
ponents in Zeeman associations have similar positional gradi- 
ents. It is possible that the very large gradient of spot LI re- 
sults from the spatial blending of two spots separated by less 
than the synthesized beamwidth. Blending of nearby maser 
spots can make it difficult to determine positional gradients 
when masers are not clearly isolated; for instance, the posi- 
tional gradient of spot n is roughly aligned with the ri dge of 
emission connecting it with spot m. iGreen etafl (120071) note 
a positional gradient aligned along a northwest-southeast di- 
rection in three of five 6035 MHz Zeeman pairs, although it 
is not clear from their text whether the ridge of spots Mm and 
Nn (which would have appeared as a single spot at their reso- 
lution) is one of them. 

4. DISCUSSION 

4.1. Zeeman Pairs 

The LCP and RCP components of 6 GHz Zeeman pairs are 
found to have a very small spatial separation. The separation 
between the centroid positions of the peak channels of the 
LCP and RCP spots of a 6.0 GHz Zeeman pair does not ex- 
ceed 1 mas except for pair Mm, the weaker component of the 
aforementioned ridge of emission in the southeast. A Zeeman 
component is identified for every 6.0 GHz maser whose emis- 
sion is not spatially blended with that of another, brighter spot. 
The 2 spots labelled Hh may constitute a Zeeman pair, but the 
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FIG. 3. — Integrated Stokes I emission and polarization of spots Mm and Nn 
in the velocity range 14.31 to 14.65 kms~'. Contours are shown as multiples 
of 4 times the rms noise. Lines indicate the direction of linear polarization. 



tion is clearly detected and found to be aligned with the ridge 
(Fig. El). 

Positional gradients are computed via weighted least- 
squares fitting of a straight line to a maser spot's centroid 
right ascension and declination in each frequency channel in 
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FIG. 4. — Relative locations of 6030 /6035 MHz masers and ground-state 
masers (in grey, fro m Fis h & Reid 2007). Adopted alignment of the two data 
sets is discussed in § 14.21 

weakness of the LCP feature prohibits measuring its velocity 
accurately enough to identify the magnetic field strength and 
direction. 

Only one Zeeman pair was identified at 6030 MHz. It 
was found to be spatially coincident with another Zeeman 
pair at 6035 MHz to within a few tenths of a milliarcsec- 
ond. The magnetic field and unshifted velocity of the 6030 
and 6035 MHz pairs agree to within 0.1 mG and 0.01 kms" 1 , 
respectively. Our m agnetic field strengths are broadly con- 
sistent with those of Desmurs & Baudry] (119981) . Differences 
on the order of 1 to 2 m G can be explained by th e lower 
spectral resolution of the Desmurs & Baudry (1998) obser- 
vations (0.11 kms" 1 channel separation =1.4 mG at 6030 
and 2.0 mG at 6035 MHz). Our magnetic field measure- 
m ents are highly con sistent with the MERLIN observations 
of iGreen et alJ ([2007), with differences well under 1 mG (no 
greater than 0.2 mG in 5 of 6 Zeeman pairs). 

4.2. Multitransition Maser Associations 

We find that nearly every 6.0 GHz maser detected in ON 1 
is found in close association with (i.e., less than 10 mas 
from) at least one ground-state maser spot (Figure |4| when 
the brightest 6035 MHz sp ot is aligned with the brightest 
1665 MHz spot detected bvlFish & Reidl d2007l) . This result 
differs from the results of IGreen et al. (2007), who find from 
lower-resolution MERLIN data that the closest 6.0/1.6 GHz 
separa tion involves a peak separation of 28 mas. iFish & Reidl 
(2007) did not obtain absolute positions for the 1665 and 
1667 MHz masers; nevertheless the separations between the 
6.0 GHz and nearest 1.6 GHz masers listed in Table|2]are al- 
ways within a few milliarcseconds, and at least one associated 
ground-state maser can be found for almost every 6.0 GHz 
maser found in this work. 

The close spatial associations allow for multitransition Zee- 
man associations to be identified. The Zeeman effect shifts 
LCP and RCP components in frequency by the magnetic field 
strength multiplied by a constant splitting factor that is dif- 



ferent for each transition. Under the assumption that mul- 
titransition associations are produced by the same cloudlets 
of material, the magnetic field strength and systemic velocity 
can be derived so long as any two Zeeman components can be 
identified from amongst all transitions seen in spatial overlap. 
Table |2] lists the magnetic fields derived from Zeeman asso- 
ciations including both 6.0 GHz and 1.6 GHz masers. There 
are four Zeeman associations in which a complete Zeeman 
pair is found at both 6035 MHz and 1665 MHz (and one ad- 
ditional pair at 1667 MHz). The magnetic fields obtained in 
these transitions are consistent with each other in all cases 
to 0.6 mG or better. The excellent agreement between mag- 
netic field strengths obtained from 1 .6 GHz Zeeman pairs and 
6.0 GHz Zeeman pairs provides strong evidence supporting 
the alignmen t between the refere nce frames of these observa- 
tions and the lFish & Reid (2007) VLBA observations. 

Furthermore, Zeeman associations consisting of a pair of 
LCP and RCP components at 6030 MHz and a single unpaired 
component at 1665 MHz show that the implied magnetic field 
at 1665 MHz (assuming that the systemic velocity is given 
by the average of the 6035 MHz LCP and RCP velocities) is 
consistent with the 6035 MHz value. In the six such associ- 
ations, the maximum difference between the 6035 MHz and 
implied 1665 MHz magnetic fields is 2.2 mG, with four asso- 
ciations showing a difference of less than 1 mG. We refer to 
these instances as "multitransition Zeeman associations" and 
argue that the above results support the validity of using these 
to derive magnetic fields. 

Of note is a potential +1.7 mG magnetic field implied by 
the overlap of a 6035 MHz RCP spot (H) with a 1665 MHz 
RCP spot at a different velocity. There is also a very weak 
6035 MHz LCP spot (h) at this position whose velocity cannot 
be determined to sufficient accuracy to identify the sign of the 
magnetic field from the 6035 MHz Zeeman pair alone. 

There does not appear to be a significant correlation be- 
tween the brighter polarization at 6035 MHz and that in a 
nearby ground-state transition, usually 1665 MHz. When 
the flux ratio between Zeeman components at 6035 MHz is 
close to unity (< 2.5), there is no correlation at all. Only two 
6035 MHz Zeeman pairs have a component flux ratio greater 
than 2.5: Hh and Kk. In both instances, the brighter Zee- 
man component at 6035 MHz is RCP, and the R CP feature at 
1665 MHz is detected, but the LCP feature is not dFish & Reidl 
120071) . In the one case where a 6030 and 6035 MHz spot 
overlap, the RCP component is brighter than the LCP compo- 
nent in both transitions, and the RCP/LCP ratio is greater for 
the 6030 MHz transition. (However, emission at 1665 MHz 
which is spatially coincident with these Zeeman pairs is seen 
only in LCP, not RCP.) There is also one instance where the 
1665, 1667, and 6035 MHz masers overlap, and the Zeeman 
component flux ratios decrease in that order. This is consis- 
tent with the fl ux ratios varying mo notonically with Zeeman 
splitting ratio (Wright et al. 2004b). There is also a similar 
trend for 6.0 GHz Zeeman pairs to have smaller separation 
between LCP and RCP components than is seen for coinci- 
de nt ground-state transitions (an effect also noted in W3(OH) 
by Fish & Siouwerman 2007) in W3(OH). These effects may 
be related if, for instance, excitation conditions (such as ve- 
locity coherence) change slightly over the physical extent of 
a Zeeman association. 

The large number of Zeeman pairs identified in transitions 
of OH allows a systemic maser velocity map of ON 1 to be 
constructed, as shown in Figure [5] Hydroxyl maser veloci- 
ties have been corrected for Zeeman splitting, which can be 
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TABLE 2 

MULTITRANSITION ZEEMAN ASSOCIATIONS 



6.0 GHz 6.0 GHz 6.0 GHz 6.0 GHz 1.6GHzLCP 1.6GHzLCP 1.6GHzLCP 1.6GHzRCP 1.6GHzRCP 1.6GHzRCP 1.6 GHz 
Zeeman B Separation v lsr R.A. Decl. vlsr R.A. Decl. vlsr B 

Group (mG) (mas) (kms^'V 1 Off. (mas) Off. (mas) (kms" 1 ) Off. (mas) Off. (mas) (kms" 1 ) (mG) b 



A 


-4.5 


0.77 


12.48 


-884.60 


707.53 


13.75 








-4.3 


B 


-4.2 


0.15 


13.94 


-855.70 


697.28 


15.15 








-4.1 


Z c 


-4.1 


0.32 


13.95 
















D 


-4.6 


0.92 


15.29 


-258.85 


100.60 


16.73 


-259.45 


100.83 


13.92 


-4.8 


E 


-4.6 


0.19 


15.32 


-249.24 


102.42 


16.73 


-249.72 


103.36 


13.68 


-5.2 


F 


-5.3 


0.87 


14.13 


-233.76 


97.58 


16.35 








-7.5 


G 


-5.1 


0.28 


14.92 


-201.74 


85.70 


15.88 








-3.2 


H 




0.92 d 


5.5 








-128.22 


994.63 


5.95 


+1.7 


I 


-12.1 


0.10 


0.14 


-50.68 


1012.23 


3.90 








-12.8 


K 


-5.2 


0.54 


14.64 








-0.02 


-0.02 


12.82 


-6.1 e 


L 


-1.7 


0.31 


13.72 


40.34 


526.99 


13.88 


39.19 


527.41 


13.21 


-1.1 


N f 


-1.3 


0.45 


14.47 


154.58 


-6.23 


15.05 


152.31 


-4.82 


14.02 


-1.7S 










161.92 h 


-13.37 


14.64 


160.688 


-13.49 


14.29 


-1.0 



NOTE. — The columns 6.0 GHz and 1.6 GHz refer to Ihe 6035 MHz and 1665 MHz transitions, respectively, except as nolcd in the table. The 1.6 GHz mascr information is laken from Tabic 2 in Fish & Reid 20071. 



Systemic velocity (corrected for Zeeman splitting). 
^ When two 1 .6 GHz spots are given, magnetic field is as derived from 1 .6 GHz Zeeman pair. When only one 1 .6 GHz spot is given, magnetic field is inferred from mttltitransition overlap. 
C 6030 MHz 

^ Interpreting the RCP and LCP features as Zeeman components. 

e Paired with a distant LCP feature to give a -2.5 mG field in Fish & Reid 2007 1. 

^ Spatially blended with Zeeman group M, which implies an almost identical magnetic field strength and systemic velocity. 
g Alternatively, the 6035 MHz spots may be grouped with (162.13.-7.20) mas LCP 1665 MHz to give a -1.0 mG field. 
h 1667 MHz 
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P 



0.5 



-0.5 



4.2. 1 . Alignment of Maser Positions 



-IS. 4 Demagnetized Velocity (km s ) 15.4- 




14 i! !« 
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6035 □ 
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Right Ascension Offset (arcscc) 

FIG. 5. — LSR velocities corrected for Zeeman splitting. Color indicates 
the central (unshifted) velocity for Zeeman pairs in each transition (also in- 
dicated with numbers), with northern masers (i'lsr "C 12.4 kms~') in black. 
Multitransitio n Zeeman ass ociations are not shown. Methanol maser data is 
taken from lRygl et alj <2010h . 



large (several km s 1 between LCP and RCP com ponents) in 
the gro und state. Me t hano l maser positions from lRygl et al.l 
(2010) and |XuetaL| (2009) are included. No such correc- 
tion is required for methanol, which is a nonparamagnetic 
molecule. There is a trend for LSR velocities to increase to- 
ward the west along the southern line of maser spots, although 
not all Zeeman pairs and methanol masers conform to this pat- 
tern. 



In aligning the 6.0 GHz masers with the lFish & Reid! (120071) 
1.6 GHz masers, we note that there is no particular rea- 
son to believe that the position of the brightest 6035 MHz 
will be precisely coincident with the location of the brightest 
1665 MHz maser to submilliarcsecond accuracy. The appar- 
ent centroids of Zeeman components may appear at slightly 
different locations. For instance, the separation between the 
components of 6.0 GHz Zeeman pairs listed in Table|2]is typ- 
ically a few tenths of a milliarcsecond. This effect is even 
larger at 1665 and 1667 MHz, where the separation between 
components of a Zeeman pair can be a few milliarcseconds. 
Assuming that there is no large-scale organization that sys- 
temically shifts LCP spots in one direction and RCP spots in 
the opposite direction, we can define a best statistical align- 
ment of the 6.0 and 1.6 GHz maser frames by finding the 
single vector (right ascension and declination) offset to add 
to the positions of one of the frames to minimize the sum 
of the squares of the 6.0 GHz-1.6 GHz position differences 
among the Zeeman associations. We use the data in Table |2j 
using the average of RCP and LCP positions within one fre- 
quency as the central location of the maser emission at that 
frequency in each Zeeman association. When only one sense 
of circular polarization is detected in a Zeeman association in 
a frequency, we use that position as the assumed central lo- 
cation of the maser emission at that frequency. We find that 
the two frames are best aligned w hen the 1665 MHz refer- 
ence feature of Fish & Reid (2007) is assumed to be 1.7 mas 
north of the 6035 MHz reference feature. In this frame, the 
median 1.6/6.0 GHz offset in the Zeeman associations is ap- 
proximately 2.2 mas. 

The offsets between 1.6 and 6.0 GHz Zeeman associations 
are consiste nt in magnitude wit h proper motions of masers in 
ON 1. The iFish & Reidl d2007l) data were observed 3.74 yr 
prior to the epoch of observations of the 6.0 GHz masers re- 
ported here. A 1.7 mas apparent offset between the epochs 
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could corres pond to a motion of 5.5 kms at a distance 
of 2.57 kpc dRvgl et allHoToh . This is comparable to the 

1.6 GHz maser motions seen in iFish & Reidl (|2007). Al- 
though the lack of phase-referencing in the ground-state ob- 
servations precludes rigorous analysis of the offsets in terms 
of proper motions, the general pa ttern is confirmed by the 

6.7 GHz methanol observations of Rvgl et al. (2010), which 
show that the relative speed of the northern and southern 
methanol masers in the plane of the sky is 9.4 km s" 1 . 

4.3. Comparison with W3(OH) 

It is instructive to compare ON 1 with W3(OH), the 
other massive star-forming region hosting a large number 
of 6.0 GHz masers observed with a modern VLBI array. 
There are numerous similarities between the properties of 
6.0 GHz masers in ON 1 and W3(OH), some of which have 
already been commented on in §§ [3] and [4] For instance, 
6030 MHz masers in ON 1 are rare and are only found 
in the presence of 6035 MHz masers. This is consistent 
with observations of W3(OH), which found that 6030 MHz 
masers are nearly always found in close spatial association 
with 6035 MHz masers an d imply similar physical conditions 
(Fish & Sjouwerman 2007). 

However, one important contrast between the OH masers in 
ON 1 and W3(OH) is that in the former, all 6035 MHz masers 
are found in close association with 1665 MHz masers. This 
stands in stark contrast with W3(OH), where many 6035 MHz 
maser spots are found without accompanying ground-state 
emission, including some (particularly in the NE and SE re- 
gions) where no ground-state masers are located within hun- 
dreds of AU (Fish & Siouwermanll2007l) . Theoretical mod- 
elling suggests that co-propagation of 1665 and 6035 MHz 
masers (as observed throughout ON 1) can occur under a 
wide range of conditions, but that the presence of 6035 MHz 
masers without 1665 MHz masers may require a highe r den- 
sity and/or higher OH column density (e.g., Gray et aflll992t 
Cragg et al. 2002). It is therefore probable that the OH masers 
in ON 1 are on the whole tracing slightly less dense material 
than in W3(OH). 

Collisions may als o play a partial role in exciting the 
6035 MHz transition (iGray et all 119921) . It is possible, for 
instance, that the masers just south of the H II region trace 
a shock front. The ap pearance of ground-state, s atellite-line 
masers in this region (Nammahac hak et "ail 120061) may also 
be indi cative of a shock dGrav et al.lll992tlPavlakis & K vlafis 
fl996aHbh . 

Methanol maser excitat ion is predicted to occur over a sim- 
ilar range of densities (Cragg et al. 2002). The tentative de- 
tection of a magneti c field strength of —18 mG in methanol 
Zeeman splitting bv lGreen et al.l d2007) may be indicative of 
a slightly higher density locally at the methanol maser site 
in the southern line of masers in ON 1. Methanol emis- 
sion is also seen intermixed with OH masers in W3(OH) 
(iHarvey-Smith & Cohenll2006l) . 

In W3(OH), 31% of the 6.0 GHz masers are 6030 MHz 
masers, which are typically weaker th an their 6035 MHz 
maser counterparts by a factor of several (Fish & Siou wermanl 
120071) . In ON 1, only 10% of the 6.0 GHz masers are in the 
6030 MHz transition. This may not be a significant differ- 
ence, due both to the smaller number of total maser spots in 
ON 1 and the fact that the 6.0 GHz masers are weaker than 
in W3(OH), meaning that a larger fraction of the presumed 
6030 MHz maser population in ON 1 ma y be below our de- 
tectability threshold. The modelling of Pavlakis & Kvlafis 



(2000) and Cragg et al. (2002) show that the parameter spaces 
under which 6030 and 6035 MHz masers are likely to be ex- 
cited in star-forming regions are very similar. For identical 
excitation conditions (as would occur if 6030 and 6035 M Hz 
masers were found in spatial overlap), Cragg et al. (2002) find 
that the brightness temperature of the 6030 MHz is typically 
a factor of a few smaller than that of the 6035 MHz maser. 
However, 6030 MHz masers are found only in the west, de- 
spite the fact that all maser clusters in ON 1 (northern, central, 
western, and the line of masers in the south) host 6035 MHz 
masers that are at least a factor of 10 brighter than our detec- 
tion threshold. It is not clear what, if any, are the differences 
among the physical conditions of these clusters that allows 
6030 MHz maser formation only in the west. 

4.4. The Structure of ON 1 

The view of ON 1 has evolved w i th inc reasingly better 
observational data. Israel & Wootten (1983) suggested that 
ON 1 was an example of isolated massive ( spectral type 
B0.5) star formati on. More recent infrared dKumar et al.1 
120021 20031 120041). centim eter, and submillimeter observa- 
tions ( Su et al.ll2004l 12009b indicate that the environment of 
ON 1 contains numerous dust cores associated with other 
intermediate-mass or high-mass star formation. Multiple 
molecular outflows and high-velocity water maser structures 
add further complexity to the surroundings of the H II region, 
suggesting that the source(s) in the H II region are gravitation- 
ally bo und to other nascent and possibly evo lved stars in the 
region dNagavama et al.l 120081; ISu et al.l 120091) . The compli- 
cated environment of ON 1 renders it difficult to interpret the 
OH masers surrounding the H II region. 

Nevertheless, based on the OH maser distribution in mul- 
tiple transitions, it is likely that the H II region is asso- 
ciated with at least two distinct excitation sources. The 
highly excited 13441 MHz transition of OH produces masers 
near both an d 14-15 kms" 1 dBaudrv & Desmursl 12001 
IFish et ali r2005). Archival VLB A observations (experiment 
code BB137) of this transition were unsuccessful in detecting 
the emission, so no interferometric map of the 13441 MHz 
masers exists. However, VLBI maps of 6.7 GHz methanol 
masers, which show a similar double-peaked spectrum at 
and 14-15 km s -1 , are consistent with emission in the northern 
and south ern groups dSugivama et al.ll20()l iRvgl et alj|20~iol: 
see also lGreen et alj|2007l) . as shown in Figure [3] The exis- 
tence of a 12 mG magnetic field, large for OH masers, in the 
north argues that the northern masers (Hh, Ii, and J) are lo- 
cated close to one of these excitation sources. The magnetic 
field measurements in the south are less conclusive, although 
if the -18 mG magnetic fiel d inferred from methanol Zeeman 
splitting dGreen et al.ll2007l) is correct, the bright 6035 MHz 
reference feature (K) may be located closest to the southern 
excitation source. We note the possibility, however, that the 
large magnetic field strengths and existence of highly-excited 
OH masers in the north and south may simply be due to dense 
molecular clumps on opposite sides of an H II region contain- 
ing only one excitation source. 

The southern masers, which constitute the bulk of the ob- 
served OH masers i n ON 1, are probably b est understood 
in the context of the Elitzur & de Jong ( 1978) model of OH 
masers forming in the compressed region between the ioniza- 
tion and shock fronts around an expanding H II region. These 
masers appear at the edge of the H II region as observed at 
A = 1.3 cm and are located not too far from the edge of the 



VLBI Observations of 6 GHz OH Masers in ON 1 



7 



dust cavity evacuated by the H II region (ISu et alJ l2009). The 
systemic velocity of the molecular cloud is fairly well es- 
tablished to be vlsr = 11-12 kms" 1 from a large variety of 
dense molecular tracers including ammonia, formaldehyde, 
CS, SiO, thermal methanol, and Class I methanol masers (e.g., 
Haschick & Baanll989t]Haschick & Hoi 1 9901: lHaschick et al. 
1990tlPlume et al.lll992ilAnglada et al.l(l9%flMacLeod et al. 
1991 iKalenskii et al.ll200Tl lArava et al.ll2002h . The line-of- 
sight velocity of the southern masers is close to the assumed 
systemic dust velocity of 1 1-12 km s" 1 but slightly redshifted 
(Figure 0, with the offset from systemic ranging from less 
than 1 to over 3 kms" 1 . This (line-of-sight) velocity offset 
is consistent in magnitud e with the 3.6 km s" 1 (transverse) 
proper motion inferred by iRygl et al.l ((2010) from methanol 
maser observations and the expansion of the H II region. Un- 
fortunately, the l inear polarization in this region (from t his 
work as well as iFish et all 120051 and iGreen et all 12007b is 
complicated enough to prevent definitive identification of the 
large-scale magnetic field direction in the plane of the sky 
along the line of southern masers. 

The northern masers are associated with the largest mag- 
netic field strength seen in the source (-12 mG at spots Ii). 
More so than in other maser clumps, the northern masers di- 
vide neatly into regions where only the LCP component is 
seen at 1665 MHz and regions where the RCP component is 
seen. After the reference feature, the 6035 MHz Zeeman as- 
sociation with the largest component flux ratio is located in 
the northern cluster. These features may be exp lainable by 
correlated magnetic fields and velocity gradi ents (Coo kTT966l) 
or by Zeeman overlap (De guchi & Watsonl 1986). There is 
also a relative dearth of linear polarization in th is region, at 
both 6035 MHz and 1665 MHz (IFish et al.ll2005l) . Indeed, the 
Zeeman splitting from the large magnetic field at spot Ii is 
larger than the maser line width, which is predicted by Zee- 
man overl ap to produce masers that are very highly circularly 
polarized (De guchi & W atson 198g). 

Understanding the northern masers may require under- 
standing both the multiplicity of massive stars in the H II re- 
gion and the velocity of the H II region. IZheng et al.l (fl985h 
measured the velocity to be +5 kms" 1 at 14.7 GHz, although 
there is reason to believe that this velocity may be signif- 
icantly blueshifted from the systemic velocity of the mas- 
sive star(s) in the H II region due to opt ical depth effects i n 
the expanding plasma (see discussion in lFish & Reidll2007h . 
w hich are v isible in the A = 3.6 and 1.3 cm continuu m maps 
of ISu et al.l (120091) . Three possibilities arise: the IZheng et al.l 
(11985|) velocity may indeed trace the velocity of a single mas- 
sive star within the H II region; there may be only one excita- 
tion source in the H II region, but its velocity is not given by 
the 14.7 GHz recombination line; or there may be more than 
one excitation source within the H II region. 

In the first case, if there is only one source of excitation 
at vlsr ~ +5 kms" 1 , the northern and southern masers are 
presumably both associa t ed wit h it. In this case, the inter- 
pretation of Fis h & Reidl (120071) is probably correct: maser 
motions are dominated by the expansion of the ionized region 
and an additional rotational component, each having a mag- 
nitude of several kms" 1 . As is usual in massive star-forming 
regions, not all OH maser velocities fit nicely into the large- 
scale velocity pattern, but it is plausible that the differences 
can be explained by projection effects and inhomogeneities in 
the density profile of the molecular material. 

In the second case, the systemic velocity of the ON 1 sys- 



tem, given by the presumed lone massive star inside the H II 
region, is likely close to the systemic molecular velocity of 
11-12 km s" 1 . Recombination line studies of other ultracom- 
pact H II regions show that centimeter-wavelength recombi- 
nation lines are typically blueshifted compared to millimeter- 
w avelength recombin ation lines (see discussion in section 3.3 
of iSewilo et al.ll2008l) . Since the optical depth of the H II re- 
gion decreases with frequency, higher-frequency recombina- 
tion lines are sensitive to the center of the H II region (where 
the massive star is likely located), while lower-frequency lines 
trace the expanding o uter portions of the H II region and are 
therefore blueshifted dBenilis & Ershovll 19831: IWelch & Marrl 
119871; iKeto et al.lll995l) . If the same effect operates in ON 1, 
the velocities of the northern masers (mostly near kms" 1 ) 
differ by more than 10 km s" 1 from the rest of the OH masers, 
the systemic velocity of the molecular material in the vicinity, 
and presumably the velocity of the star in the center of the 
H II region. It is possible, then, that the northern masers trace 
an outflow, possibly with a helical magnetic field. Ground- 
state RCP and LCP features have a tendency to segregate on 
opposite sides of the presumed structure, which could be ex- 
plained if the magnetic field wraps around the outflow. Unfor- 
tunately the magnetic field direction and strength at spot Hh 
are poorly determined due to the low signal-to-noise ratio of 
the LCP spot, so it is not possible to confirm that the polarity 
of the magnetic field reverses across the northern structure. 

In the third case, there are at least two massiv e stars in the 
H II re gion, each with associated OH masers. The Zheng et al. 
( 1985) recombination line velocity may represent an average 
of the velocities of the two sources, although optical depth 
effects may still be important. Based on both the recombina- 
tion line velocity and the fact that the northern masers (un- 
like the southern masers) are projected atop the H II region, 
it is probable that the northern excitation source is in front of 
the southern excitation source. Either the expansion/rotation 
or outflow interpretation of the maser distributions may ap- 
ply, depending on the relative velocities of the two excitation 
sources. 

In contrast with observations of W3(OH), the 6.0 GHz 
masers in ON 1 do not clearly trace additional structures that 
were not visible at 1.6 GHz, limiting their ability to provide 
additional information on the context of the various maser 
clusters in the source. Further understanding of how the maser 
clusters are associated with one another will require higher- 
resolution observations of the nonmasing material. Sensitive 
observations of high-frequency recombination lines in the H II 
region may be illuminating regarding the internal structure of 
the ionized region, especially if done with enough spatial res- 
olution to understand the internal dynamics of the ionized re- 
gion. ALMA observations of dense molecular tracers may 
also assist in providing the context in which star formation is 
taking place in ON 1 . 

5. CONCLUSIONS 

We have imaged the 6030 and 6035 MHz OH masers in 
ON 1 at VLBI resolution. The distribution of these excited- 
state masers is similar to that of the ground-state masers. 
Unlike in the similar massive star-forming region W3(OH), 
6.0 GHz masers are not found in the absence of 1665 MHz 
masers, perhaps suggesting that ON 1 does not have ana- 
logues of the highest-density knots found in W3(OH). The 
6.0 GHz masers are spatially coincident with 1665 MHz 
masers to within a few milliarcseconds. Magnetic fields 
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strengths obtained from Zeeman splitting at 6.0 and 1.6 GHz 
are usually consistent to better than 1 mG, suggesting that 
multitransition Zeeman associations are acceptable for obtain- 
ing estimates of the local magnetic field. 

Our observations confirm the existence of a strong 
(-12 mG) magnetic field among the northern, blueshifted 
masers, as s uspected from early EVLA observations of ON 1 
dFishl 12007). The large magnetic field here, as well as the 
distribution of methanol maser spots and presumed distribu- 
tion of the highly-excited 13441 MHz OH masers, may indi- 
cate that the northern and southern masers are excited by two 
different sources within the H II region or that the northern 
masers trace and outflow. However, other scenarios, such as 



large-scale expansion and rotation of all groups of masers in 
ON 1, cannot be ruled out from the data at hand. The overall 
structure of ON 1 remains uncertain in the absence of sensi- 
tive, high-resolution observations of the nonmasing material. 

The National Radio Astronomy Observatory is a facility 
of the National Science Foundation operated under coopera- 
tive agreement by Associated Universities, Inc. The European 
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